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FOREWORD 


This  report  describes  the  results  to  date  of  an  analysis  conducted 
by  the  University  of  Tennessee  Space  Institute  under  U.S.  Army  Contract 
No.  DAAHOI-69-C-1357.  The  contract  was  initiated  under  DA  Project  No. 
1M2326XXA206  and  AMC  Management  Structure  Code  No.  522C  .11 .14800.  The 
technical  effort  was  performed  between  March  1970  and  April  1971  under 
the  direction  of  the  Aerodynamics  Group,  Aeroballistics  Directorate, 

U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama.  The  Army  technical 
representative  was  Mr.  D.  J.  Spring. 
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Section  I.  INTRODUCTION 


The  progress  made  since  publication  of  our  last  summury  report  |1) 
is  summarized  and  reported  herein.  Considerable  knowledge  was  obtained 
in  the  inviscid  flow  field,  which  is  necessary  before  one  can  study  in 
detail  the  viscous  flow  field  and  interaction  on  various  bodies  of 
revolution  including  small  angles  of  attack.  The  studies  on  viscous 
effects  are  preliminary  ones;  however,  some  interesting  results  as  well 
as  the  direction  of  future  investigations  ar;  pointed  out.  Since  this 
is  a  summary  report,  most  details  were  omitted;  however,  a  series  of 
publications  which  present  the  calculation  procedures  as  well  as  detail 
discussions  are  referenced. 

The  major  ideas  and  typical  results  on  the  inviscid  and  the  viscous 
flow  studies  are  presented  in  Sections  II  and  III,  respectively.  The 
analysis  on  inviscid  transonic  flow  field  around  an  ogive-cylindrical 
body  at  zero  angle  of  attack  (2]  has  been  extended  to  include  the  effect 
of  cross-flow  caused  by  small  angles  of  attack.  The  results  are  very 
encouraging  and  check  reasonable  well  with  experiments  including  recent 
investigations  by  AMICOM.  The  flow  over  a  conical  nose  shape  has  been 
solved  using  the  integral  approach  for  the  entire  transonic  flow  regime. 

It  will  be  shown  that  this  is  a  higher  order  approximation  method.  The 
analysis  can  take  into  account  the  cone  at  small  angles  of  attack  with 
good  agreement  with  the  available  data.  This  study  has  been  extended 
to  calculate  the  flow  field  along  a  cone-cylindrical  body  at  zero  angle 
of  attack  where  the  method  used  for  the  cylindrical  body  portion  has 
been  reported  in  previous  publications  [  1 ,  3]  . 

The  effect  of  viscous  and  inviscid  flow  interaction  is  studied  by 
considering  the  boundary  layer  build  up,  including  transition  and  through 
a  normal  shock  wave,  and  its  results  upon  the  inviscid  flow  field.  (The 
interaction  between  the  potential  solution  and  the  boundary  layer  solution 
is  studied  through  the  iteration  process.)  Remarks  concerning  the 
usually  accepted  "equivalent  body"  concept  are  noteworthy ,  especially, 
if  separation  is  involved.  A  preliminary  study  on  the  separation  caused 
by  the  exhaust  jet  plume  indicates  that  the  transonic  separation  could 
be  very  different  from  that  of  the  lower  supersonic  case.  Moreover,  if 
the  separated  region  is  small,  the  usually  accepted  "equilibrated"  free 
shear  layer  analysis  is  not  applicable. 

A  very  basic  study  program  on  the  transonic  separation  problem  has 
been  initiated.  Studies  involving  problems  in  improving  the  inviscid 
flow  solutions  are  continuously  in  progress. 
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Section  II.  INVISCID  FLOW  FIELD  AROUND  VARIOUS 
BODIES  OF  REVOLUTION 


The  nonlinear  correction  theory  with  stretching  procedure  has  been 
extended  to  a  lifting  ogive- cylindrical  body.  The  solutions  were 
obtained  over  the  entire  transonic  flow  regime.  The  major  ideas  and 
some  results  are  presented  in  Paragraph  1.  The  integral  equation 
approach  to  the  flow  field  along  a  conical  nose  is  discussed  in 
Paragraph  2.  The  complete  solution  over  a  nonlifting  cone-cylindrical 
body  is  discussed  in  Paragraph  3. 


1.  Flow  Field  Around  Ogive-Cylindrical  Bodies  at  Small  Angle  of  Attack 

a .  Basic  Consideration 

It  has  been  shown  [1,  2]  that  Hosokawa's  nonlinear  correc¬ 
tion  theory  [4,  5]  with  a  stretching  procedure  [2]  can  be  applied  to  a 
cylindrical  body  with  an  ogive  nose  at  zero  angle  of  attack.  The  surface 
pressure  distributions  as  well  as  shock  wave  locations  over  various  body 
configurations  have  been  shown  to  agree  fairly  well  with  experiments. 
Encouraged  by  these  good  agreements,  an  extension  to  include  the  body 
at  small  angle  of  attack  has  been  attempted.  The  following  discussion 
only  presents  the  major  idea  and  the  detailed  approach  is  in  preparation 
as  a  separate  publication*. 

The  appropriate  inviscid  transonic  flow  small  perturbation  equation 
can  be  written 

(  1  ‘  Ml)  ^xx  +  7  ^  (rir  )  +  ^  =  <7  +  Dm2  <i>x^xx  (l) 

in  which  the  velocity  potential,  <$,  is  normalized  with  respect  to  the 
free-stream  velocity  .  The  coordinates  and  the  geometry  are  illus¬ 
trated  in  Figure  1.  A  solution  must  satisfy  the  boundary  conditions  of 
tangential  flow  at  the  body  surface  with  the  perturbation  velocities 
vanishing  at  infinity. 

The  tangential  boundary  condition  can  be  approximated  for  a 
slender  body  by 


*Wu,  J.  M.  and  Aoyama,  K.,  Transonic  Field  Around  Ogive-Cylinder  at 
Small  Angle  of  Attack,  U.S.  Army  Missile  Command,  Redstone  Arsenal, 
Alabama,  Technical  Report  (to  be  published  soon),  1971. 
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Nomenclature  of  Tangent  Ogive  Body 


Tor  +  Sin  0!  Cos  9  1  , 

L  1  J  on  body 


dR 

dx 


(2) 


where  a  and  R  are  the  angle  of  attack  given  in  radians  and  the  normalized 
body  radius  with  respect  to  the  total  body  length  L,  respectively. 

In  Hosokawa's  nonlinear  correction  theory  [5],  the  nonlinear 
solution  4  of  Equation  (1)  is  assumed  to  be  formed  by  the  linearized 
solution  <>  of  Oswatitsch  and  Keune  [6]  type  and  a  correction  function  g. 
it  is  a  solution  of  the  equation  which  reads 


Id/  \ 

7“(r*r) 


1  d^$ 

f2  dt?2 


(3) 


where  the  contribution  from  the  azimuth  direction  is  included  because  of 
the  cross  flow  induced  by  the  small  angle  of  attack.  K  is  an 
acceleration  like  term  and  has  been  assumed  to  be  a  constant.  This 
assumption  leads  to  the  so-called  "parabolic"  method  known  to  the 
technical  community.  Hosokawa  assumed  that 


4>  =  <t>  +  g 


(4) 


where  g  is  the  correction  function.  A  further  study  on  this  g  function 
is  underway  and  will  be  in  a  separate  technical  report*. 

By  Equations  (1),  (3),  and  (4),  and  by  order  of  magnitude  arguments 
[ 5] ,  it  has  been  found  that  the  g  function  has  to  satisfy  the  following 
equation: 


+  (1-7) 


-  (1+7)  M2 


'xx 


(5) 


The  correction  function  g  may  be  understood  to  be  a  function  of  x  with 
r  and  9  as  parameters,  i.e.,  g(x;r,9).  Integration  of  Equation  (5) 
gives 


*Wu,  J.  M.,  Some  Remarks  on  Transonic  Flow  Nonlinear  Correction  Theory 
(in  preparation) . 
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gx  (x;r,g) 


(6) 


1-M 


QQ 

(1+y  )Mi 


where 


Y(x) 


l-M* 

-  _ _ 

2 

_  o 

X 

r 

"t  -  K 

(14  7)M2 

*  Cm 

y* 

X 

C'+Dm2 

*x  dx 


in  which  x*  is  the  unknown  parabolic  point  where  the  coefficient  of  £xx 
in  the  nonlinear  transonic  equation  vanishes. 

The  double  sign  of  Equation  (6)  is  determined  according  to  the  first 
term  on  the  right  side  of  Equation  (6)  [2,  5],  i.e., 


UMi  > 

(7+1  )M« 


0 


(7) 


It  is  suggested  that  the  parabolic  point  x*  is  determined  such  that 
the  x  component  of  perturbation  velocity  determined  from  the  linearized 
transonic  equation,  gives  the  nonlinear  solution  at  that  point,  i.e., 


4>x  (  x*»r»e  )  ■ 


4X  ) 


(7+1  )M  2 


(8) 


This  satisfies  the  condition  that  the  acceleration  is  continuous  in  the 
accelerated  flow  regime.  From  comparison  of  Equations  (1)  and  (3)  with 
Equation  (8),  one  obtains, 


$ xx  (**;r>o)  *  - t- 

( 1+7  )M  » 


(9) 


Therefore,  Equations  (8)  and  (9)  are  the  two  equations  which  determine 
the  two  unknowns,  K  and  x*.  The  nonlinear  solution,  can  thus,  be 
constructed  according  to  Equation  (4)  with  Equation  (6)  as  follows: 


l-M: 


- ±  AW 

(7f  1)M  J 


(10) 
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where 


Y(x)  =  -2  «x(x*)  Ux(x)  -  *x(x*)  1  +  — ^ - -  [*<x>  -  *(x*)l  * 

V  '  1  J  (7+DMj  L  J 


The  linear  perturbation  potential  is  assumed  to  be  expressible  as 
the  sum  of  two  parts  because  of  the  axial  flow  and  the  cross  flow.  This 
assumption  is  valid  only  if  the  angle  of  attack  is  small  and  no  boundary 
layer  separation  is  present  in  the  flow  field.  This  gives 


$  (x,r,0)  =  0a  (x»r)  +  fc  (x>r»e) 


(ID 


where  <t>a  and  $e  denote  the  axial  and  cross  flow,  respectively. 


By  the  linearized  small  perturbation  Equation  (3),  the  governing 
equations  for  <t>a  and  <tc  read 


(l  •  "*)  **x*+  H  (C  ‘ar)‘  11  ♦> 


(12) 


and 


(‘  ■"-)  %x  +  r^(r‘cr)  +  rT 


o2  «e 


K  9, 


(13) 


The  solution  for  the  axial  flow  <j»a  has  been  discussed  in  detail  by  Wu 
and  Aoyama  [2]  .  The  discussion  on  the  solution  for  the  cross  flow  <J>C 
is  given  by  Liepmann  and  Roshko  [7] .  The  simplified  cross  flow 
solution  reads 


♦c  (x>r»0) 


cos  a  Lh 

dr 


Sin  a  Cos  d 


r2(x) 


(14) 


The  last  expression  is  known  as  the  Munk-Jones  cross  flow  term. 

The  boundary  condition.  Equation  (2),  can  be  rewritten*  for  the 
linear  solution  as  follows: 


*For  details,  see  Wu  and  Aoyama,  loc.  cit. 
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body 


dR 

dx 


(15) 


a 


r 


and 


K 


+  Cos  9  Sin 


a)l 


body 


(16) 


Now,  if  the  angle  of  attack  is  small,  it  is  reasonable  to  assume 
that  the  linear  perturbation  velocity  may  be  expressed  as 


<t>  (x,r,0)  =  <t>0  (x,r)  +  A  ♦  (x,r,0) 


(17) 


where  the  subscript  o  denotes  the  basic  solution  at  zero  angle  of  attack, 
i.e.,  the  axial  flow  solution.  Thus,  by  comparison  of  Equations  (11) 
and  (17),  one  identifies  A  <t>  with  the  Munk-Jones  cross  flow  term. 

However,  it  should  be  noted  that  is  not  the  same  solution  of  the 

axisymmetric  equations  discussed  in  reference  [2],  because  of  the  shift 
in  the  sonic  point  location  and  the  change  in  the  constant  K. 

By  the  conditions  set  out  in  Equations  (8)  and  (9),  one  obtains 
<t>°x  (x0*  +  Ax*,  Ko  +  Ak)  +  A*x  (xo*  +  Ax*,  KQ  +  AK) 

=  (1"Mi)  /  U+7)M*  (18) 

and 

Vx  (X°*  +  Ax*’  K°  +  AK)  +  A  4,3131  (X°*  +  Ax*’  K°  +  AK) 

=  (ko  +  AK 

(1+7>mJ 


(19) 


where  xQ*  and  Kq  are  the  solutions  based  on  the  axisymmetric  solution 

[2]  . 
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The  nonlinear  perturbation  velocity  0X,  finally  reads 


1-M' 


(7+D  M2 


±  yY0  +  AY 


(20) 


where  YQ  is  due  to  the  axial  flow  and  AY  is  due  to  the  small  angle  of 
attack.  The  expressions  for  Y0  and  aY  are: 


Y0  (x) 


i-mJ 

-2  - 

(1+7  )*** 


d+7  )M  1  J 


+ 


2Ko 

(1+7  )M2 


*o<x> 


(21) 


and 

2 

ay  (x)  =  -2  L~M”  A  <t>x(x)  +  __2AK _  <to(x) 

d+7  )Mi  (1+7  )M  2 

+  2K°— -  [*A  0(x)  -  AK«0  (x0*.  K0) 

d+7  )Mi  L 

-  Ax*  (x0*,  K0  )]  (22) 


therefore,  the  correction  of  Y (x)  due  to  the  angle  of  attack  can  be 
performed . 


b .  Discussions  of  Solutions  and  Comparisons  with  Data 

Without  going  into  details  (which  will  be  published  as  a 
senarate  report*),  some  typical  solutions  are  presented  in  Figures  2 
through  7 . 

The  pressure  distribution  along  the  d  =  90  degree  surface  changes 
with  angle  of  attack  as  illustrated  in  Figure  2.  An  increase  in  angle 
of  attack  increases  the  shock  wave  strength  as  expected.  In  the 


*Wu  and  Aoyama,  loc.  cit. 
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Figure  3.  Mach  Number  Contours  on  Leeward  and  Windward 
Planes  for  5-Degree  Angle  of  Attack  at  ■  0.975 
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Figure  4.  Mach  Number  Contours  on  Leeward  and  Windward 
Planes  for  5-Degree  Angle  of  Attack  at  =  1.00 
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Figure  6.  Comparison  of  Present  Theory  with  Experiments 
Leeward  Side  lor  2-Degree  Angle  of  Attack  at  ^  ■  0.95 
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1 


axisymmetric  flow  case  [2],  the  pressure  coefficient  along  the  cylindri¬ 
cal  body  vanishes.  However,  this  is  not  so  in  the  angle  of  attack  case 
because  of  the  cross  flow. 

The  entire  flow  field,  disturbed  by  the  presence  of  a  body,  is 
illustrated  in  Figures  3  and  4  for  a  subsonic  and  a  sonic  free  stream 
case,  respectively,  for  the  plane  including  the  leeward  and  the  windward 
variations.  The  flow  pattern  on  the  leeward  side  changes  more  gradually 
compared  with  the  windward  side.  It  can  be  seen  in  Figure  3  that  the 
movement  of  the  sonic  point  is  more  sensitive  to  the  movement  of  the 
shock  location  at  angle  of  attack. 

Figure  5  presents  an  indication  of  the  circumferential  pressure 
distributions  at  various  stations.  The  suction  reaches  a  maximum  near 
the  shoulder.  The  major  lift  component  comes  from  the  front  portion  of 
the  body.  On  the  cylindrical  portion  of  the  body,  the  cross  flow  contri¬ 
bution  is  more  significant. 

A  comparison  was  made  with  recent  AMICOM  data  for  vurious  angles  of 
attack  over  the  entire  transonic  flow  region  and  the  agreement  is  very 
good.  Two  typical  examples  of  windward  and  leeward  pressure  distribu¬ 
tions  compared  with  the  theory  are  presented  in  Figures  6  and  7, 
respectively.  More  comparisons  will  be  included  in  future  publications. 
The  Inclusion  of  the  viscous  effect  will  be  discussed  in  Section  III. 


2.  Flow  Field  Along  Nom  Conos 

a .  Basic  Consideration 


The  analysis  of  the  inviscld  flow  around  a  conical  body 
at  various  transonic  speeds  has  been  performed  with  much  success.  This 
is  very  encouraging,  in  view  of  the  fact  that  the  existing  available 
techniques  are  restricted  to  a  free  stream  of  Mach  number  1  (or  a  very' 
minor  perturbation  from  1).  Moreover,  the  present  method  can  also  take 
care  of  the  small  angle  of  attack  case. 

The  appropriate  small  disturbance  transonic  flow  equation  is  Equation 
(1),  given  in  Paragraph  l.  For  the  present  case,  all  the  lengths  are 
normalized  with  respect  to  the  axial  length  of  the  cone  (Figure  8). 

Now,  in  case  of  the  transonic  flow  past  a  cone,  it  is  a  well  known  fact 
that  the  sonic  velocity  is  fixed  at  the  shoulder,  and  the  flow  is  sub¬ 
sonic  over  the  entire  conical  forebody,  as  long  as  the  angle  of  attack 
is  small  and  the  shock  wave  is  not  attached  to  the  body. 


Rewriting  the  nonlinear  transonic  equation,  Equation  (1),  as 


(l-M2J  -  (l  +  y)Mi 


U  (1  * ) 
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then,  the  above  observation  for  subsonic  flow  over  the  entire  forebody 
suggests  that  the  coefficient  of  4*xx  has  to  be  positive  and  nonzero. 

This  inspired  the  following  linearization  technique: 


(l-M2w)  -  (1+7)hI  *x  =  f2(x)  (23) 

where  f(x)  is  an  unknown  function  which  will  be  determined  later.  This 
approach  is  somewhat  similar  to  that  devised  by  Cole  and  Royce  | 8]  in 

which  •:/  x  is  approximated  but  "2:/"X2  is  not. 

Now,  by  introducing  a  new  independent  variable  :  such  that 


f 


dx 

d‘ 


(24; 


the  original  equation,  Equation  (1*;  is  reduced  to 


■*e 


t 


2 


(23; 


which  is  a  Laplace  equation,  and  ;  is  the  linearized  velocity  potential 
normalized  with  respect  to  the  free  stream  velocity. 

As  approximated  In  Paragraph  1,  assume  that  the  perturbation 
velocity  potential  :  is  composed  of  the  axial  flow  contribution  :a  and 

the  cross  flow  term  jc  is  due  to  the  small  angle  of  attack.  This 
gives 


i  ('  ,r,  ■)  -  :a  (*  ,r)  +  :c  (*  ,r , 


(26) 


The  appropriate  boundary  conditions  for  :a  and  :c  are 


<c  (Sri  )  ■  Sin  Cos  R2(f)/r 


(2?) 


and 


r) 


body  surface 


(28) 
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where  ;  is  the  half-angle  of  the  cone  measured 
The  fundamental  solution  to  Equation  (25) 


in  radians  . 
for  ?a  is 


<F,r) 


1 

4  n 


i*  F(t) 

f  ^(f-t)2  +  r2 


dt 


(29) 


where  ?*  is  the  sonic  velocity  point  in  the  transformed  plane  corresponding 
to  the  shoulder  location,  and  F(t)  is  the  source  distribution  along  the 
body  axis  determined  by  the  boundary  condition  of  the  problem, 
respectively  . 

Without  going  into  details  (the  details  will  be  published  as  a 
separate  report*),  it  is  found  that 


o  <t>a 


or 


F(|) 

2  rt  r 


(30) 


where  the  source  distribution  F(£)  satisfies 

F(fj)  =  2n  R(fJ  =  2jt&2  X (| )  . 

The  relation  between  X  and  ?  is,  from  Equation  (24), 


X(F) 


f(t)  dt 


Therefore , 


F(‘) 


f(t)  dt 


(31) 


(32) 


(33) 


and,  with  Equations  (29)  and  (26),  one  obtains  the  velocity  potential 
as , 


*Wu,  J.  M.,  and  Aoyama,  K.,  On  Transonic  Flow  Past  a  Cone-Cylinder  Body 
with  and  Without  Angle  of  Attack,  U.S.  Army  Missile  Command  Technical 
Report,  Redstone  Arsenal  (to  be  published  soon),  1971. 


lb 


■  -4/ 


*  f 


+  Sin  Q  Cos 


*  ^ 


f (s)ds 

f-  \  2  I  y.2 


(B-t)*  +  r 


f(t)  dt 


To  obtain  the  expression  for  f(f)>  following  Equation  (34),  one  obtains 


'(f-f*)2  +  r2 


2  2 

(5-t )  +  r 


+  2  —  Sin  a  Cos  f  (e ) 

r 


f(t)  dt 


However, 


f2(') 


(■  +  l)Mi 


where  the. relation  between  •:/  x  and  f(c)  in  Equation  (23)  was  employed, 
Then,  Equations  (35)  and  (36)  arc  reduced  to  the  following  integral 
equation  along  the  body  surface, 

[  f(o]3  -  f(l-Mj)  -  2  ;  (1+ >  )  M2  Sin  ,  Cos  f(0 

+  \  b2  (1+  y)  M2  [  f  LLU  Jt.  .  1  1 

[i  yC-t)Z  +  R2(‘)  +  K 2 ( *  ) J  . 


The  function  f(‘)»  which  behaves  as  the  square  root  of  the  perturba¬ 
tion  velocity,  must  satisfy  the  boundary  conditions  at  the  apex  and  at 
the  shoulder.  At  the  apex  where  x  ■  0,  the  stagnation  condition 
requires 


Ueo  +  Ho  *x 


0 


or 


*x  =  -i 


(38’) 


From  Equation  (23)  and  by  Equation  (38'),  the  solution  should  asymptoti¬ 
cally  behave  as 


f(o) 


(39) 


At  the  sonic  point  where  |  =  <•*,  one  expects 


f(|*)  =  0 


(40) 


By  employing  Equation  (32),  the  sonic  velocity  location  is  at  the 
shoulder  and  is  determined  by  the  condition 


e* 

f  f(t)  dt  a  1 

o 


(41) 


With  these  boundary  conditions,  the  nonlinear  integral  equation, 
Equation  (37),  is  solved  for  f  (?,  )*.  The  velocity  potential  is  then 
obtained  by 


A 

X 


(?) 


(l-Hj)-  [  f(pl  2 

(y-f  1)M  2 


(42) 


Then,  x  (f)  is  related  to  f  (p)  by  the  condition  of  Equation  (32).  The 
pressure  coefficient  on  the  body  surface  is  determined  by 


Cp  -  -2  *x  -  *2  +  (  1-4  Sin2  ■> )  Sin2  a 


(43) 


*For  the  details  of  solving  Equation  (37),  see  Wu  and  Aoyama, 
loc .  c it . 
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b . 


Discussions  of  Solutions  and  Comparisons  with  Data 


The  limited  available  theories,  which  were  developed  by 
quite  different  approaches  for  sonic  flow  over  a  cone-cylindrical  body, 
are  compared  to  the  present  theory.  The  available  theories  are  the 
local  linerization  method  of  Sprieter  and  Alksne  [9],  Yashihara's 
numerical  approximation  [10],  and  Leiter  and  Oswatitsch's  time  dependent 
characteristic  method  [11].  These  are  compared  to  the  original  data  of 
Page  [12]  and  the  present  theory  in  Figure  9.  The  present  result  and 
Leiter  and  Oswatitsch's  fall  very  close  to  each  other  and  agree  very 
well  with  the  data. 

It  may  be  noted  that  the  present  theory  predicts  reasonable  values 
along  the  center  portion  of  the  cone.  The  theory  predicts  a  somewhat 
poorer  solution  near  the  shoulder  because  of  the  approximation  used  in 
solving  the  integral  equation  for  f(t)  (discussed  in  detail  in  a  report 
to  be  published).  The  poor  agreement  near  the  apex  point  is  caused  by 
violating  the  small  disturbance  assumption.  Treating  the  apex  and  the 
shoulder  as  singular  points  and  matching  of  the  asymptotic  solutions  agree 
witli  the  present  result  away  from  the  singularity. 

The  applicability  of  the  present  method  to  flows  other  than  the  sonic 
flow  are  demonstrated  in  Figures  10,  11,  and  12.  The  agreement  of  the 
data  for  the  slightly  subsonic  flow  of  Figure  11  and  slightly  supersonic 
flow  of  Figure  12  is  very  encouraging. 

Typical  results  for  bodies  at  small  angles  of  attack  are  shown  in 
Figures  13  and  14*  It  can  be  seen  that  the  pressure  change  with  azimuth 
angle  is  rather  significant  on  the  windward  side.  In  Figure  15,  the 
variation  of  pressure  distribution  resulting  from  different  angles  of 
attack  is  shown.  A  comparison  with  data  is  also  given  in  Figure  16  for 
a  slightly  blunted  cone  at  a  6-degree  angle  of  attack. 


3.  Flow  Field  Along  Cone-Cylindrical  Bodies 

The  analytical  treatment  of  the  flow  downstream  of  the  sharp 
shoulder  of  a  cone-cylindrical  body  configuration  has  been  studied  and 
reported  [3].  Therefore,  by  matching  the  cone-solution  as  described  in 
Paragraph  2  to  the  developed  local  two-dimensional  approximation  method, 
it  is  possible  to  obtain  the  complete  "on  body11  solution  for  a  cone- 
cylindrical  body  at  zero  angle  of  attack.  The  restriction  to  zero  angle 
of  attack  is  because  the  work  derived  and  reported  in  reference  [3]  was 
based  upon  the  axisymmetric  condition. 

From  the  conical-nose  portion  of  the  analysis,  the  local  Mach 
number  at  the  shoulder  has  been  taken  as  unity  in  transonic  flow.  A 
Prandtl-Meyer  expansion  around  the  corner  is  then  assumed.  The  local 
two-dimensional  approximation  method  permits  transformation  between  the 
axisymmetric  and  the  two-dimensional  bodies.  The  final  result  for  the 
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igure  9.  Comparison  of  Various  Methods  with  Data  for  7-Degree  Cone  at  Sonic  Speed 
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Figure  10.  Comparison  of  Present  Method  with  Local 
Linearization  Method  and  Data  at  Sonic  Speed 


Figure  11.  Comparison  of  Present  Method  with  Data 
for  10-Degree  Cone  at  *=  0.9 
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Figure  12.  Comparison  of  Present  Method  with  Data 
for  15-Degree  Cone  at  M  =  1.10 
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perturbation  velocity  along  the  cylindrical  portion  of  the  body  calculated 
from  the  shoulder  downstream  can  be  expressed  as  [  3 ) : 


x  -  c 


x  5 

and  C  is  given  by 

c2  .  25  (1+  r)^/3 

2  1-lg 


M 


2/3 


<l+y)M2 


(x  -  c)2 
(l+7)^/3 

i-»U2 

U+r)H5 

1/2 

i  1-M2  1 

-  oo  . 

[5/  l-M2 

2 

\ 

2  o  + 

(1+  7  )  *£  | 

4\(l+7)^/ 

2/3 

(  \  1 

/  3t> 

V3 

\2v1+7  y 

l  U2 

where  •  is  the  semiapex  angle  of  the  cone.  The  pressure  coefficient 
along  the  cylinder  reads 


CP  C2  - 


2 

5 


(x  -  c) 

v/(l+7)  M2w/3 


_1 

25 


(l+7)  M2/3 


d+7)^ 


1/2 


Some  typical  calculated  results  compared  with  experiments  are  given 
in  Figures  17  through  19  for  subsonic,  sonic,  and  supersonic  freestream 
Mach  numbers,  respectively.  The  shock  condition  and  the  possible 
boundary- layer  separation  for  subsonic  freestream  Mach  number  cases  will 
be  discussed  in  more  detail  in  a  report  to  be  published  on  cone- 
cylinders  in  transonic  flow*. 


*Wu  and  Aoyama ,  loc ,  cit . 
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Section  III.  CONCERNING  VISCOUS  EFFECTS 


This  part  of  the  report  discusses  two  main  topics.  Paragraph  1  deals 
with  the  distortion  of  the  potential  flow  field  due  to  the  development  of 
the  boundary- layer  on  the  body  surface,  while  Paragraph  2  concerns  the 
interaction  of  a  jet  exhaust  plume  with  the  body  flow  field. 

The  work  described  is  the  first  stage  in  the  determination  of  the 
entire,  viscous,  flow  over  a  missile  configuration  at  transonic  speeds. 

It  will  be  shown  that  such  calculations  can  be  performed  adequately;  with 
the  exception  of  certain  vital  components.  Areas  where  little  success 
can  be  claimed,  relate  to  the  turbulent  boundary-layer  close  to  separa¬ 
tion  in  a  subsonic  flow,  and  to  a  turbulent  boundary- layer  developing 
on  a  body  containing  discontinuities  in  slope.  Work  is  in  progress  to 
redress  these  deficiencies. 


1.  Studies  on  Viscous/lnviscid  Interactions 

Writing  the  Navier  Stokes  equations  (under  the  assumption  of 
constant  coefficient  of  viscosity) 


+  —  grad  P  e  v  [curl  curl  ^  +  4  grad  A 

Dt  P  L  J 


d 


+  div  (pq)  =  0 


(44) 


where  for  the  general  vector  A 


D  A  <3  A  ,1  j  »2 
dT  a  ST  +  2  grad  A 


A  x  curl  A 


and  A  is  the  dilatation  div  q;  we  consider  only  steady  flow,  wherein  all 
temporal  derivatives  are  zero. 


Then,  writing 


q  - 

% 

+ 

/-N 

1—1 

sj 

+  •  •  • 

p  ~ 

Po 

+ 

c  <\) 

pl 

+  .  .  . 

p  ~ 

po 

+ 

c  OO 

P1 

+  .  .  . 

(45) 
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where  \  ■  1/Re  with  Re  the  Reynolds  number,  and  substitution  into  Equation 
(44)  gives  in  the  limit  \  -*  0 


D  q0  i 

—  +  -  grad  P0  =  0 


div 


■  0 


(46) 


The  subscripts  o,  1,  and  2  indicate  the  order  of  the  solutions. 


These 


curl 


equations  are  the  classical  potential  flow  equations  -  if 
— > 

0,  so  that  q0  s  grad  q»0 


when 


0 


div  grad  -,,0)  =  0 


(46a) 


These  equations  have  been  solved  for  transonic  flow  as  discussed  in 
Section  II  and  previous  publications  [1,  2], 

Transforming  the  independent  variables 


-1/2 

x  =  X;  y  =  Y\ 
where  X,  Y  ~  0(1),  and  placing, 


p  ~  p  1 
ro 

+ 

£>1<\) 

+  .  . 

• 

i 

i’~  f'o 

+ 

B1(\)  i 1 

+  .  . 

(47) 

with  similar  expressions  for  the 

velocity  components, 

gives  in  the  limit 

(\  -»  0) 
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div 
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+ 


v 


O^u. 


(48) 


as  the  first  approximation  to  the  inner  expansion  (i.e.,  the  classical 
boundary- layer  equations).  The  solution  to  these  equations  will  be 
indicated  in  Paragraph  2  for  laminar  and  turbulent  flow. 

liquations  (46)  and  (48)  represent  the  classical  solutions  to  the 
flow  over  a  body.  However,  these  solutions  do  not  introduce  any  inter¬ 
action  between  the  boundary- layer  development  and  the  external  flow. 

Only  through  the  evaluation  of  (at  least)  the  second  terms  in  the  sequences 
[Equations  (45)  and  (47)]  can  this  interaction  be  accounted  for. 

To  some  approximation,  it  is  possible  to  evaluate  the  interaction 
by  utilizing  the  equivalent  body  concept.  Here  the  boundary- layer  dis¬ 
placement  surface  is  added  to  the  geometric  body,  thus  generating  an 
"equivalent  body".  The  potential  flow  is  then  determined  for  this 
equivalent  body  and  an  iteration  is  established  between  the  body  shape 
and  the  boundary- layer  development. 

It  should  be  noted  that  for  steady  flow, 


1 

2 


grad  q 


2 


X  w 


p 


V 


curl  curl  q 


so  that  only  if  v  =  0  and  the  vorticity  vanishes  ■  0)  is  the  pressure 
gradient  aligned  with  the  velocity  gradient.  The  assumption  of  zero 
vorticity  implies  an  irrotational  approaching  flow  and  that  any  shock 
waves  in  the  flow  are  sufficiently  weak  for  their  entropy  gradients  to 
be  ignored.  In  viscous  flow  this  is  not  so  and  the  additional  terms, 
v  curl  w,  etc  must  be  included.  In  the  equivalent  body  technique, 
these  terms  are  forced  to  be  zero  outside  the  displacement  surface  and 
the  streamlines  in  the  boundary  layer  do  not  "match"  those  of  the 
external  flow.  The  potential  flow  past  an  equivalent  body  does  not 
model  the  real  flow  [13].  However,  unless  the  flow  is  near  to  separation, 
the  correction  is  small  (d5*/dx  small)  and  the  equivalent  body  technique 
gives  a  good  result. 

Figure  20  presents  a  typical  calculated  boundary- layer  development 
for  a  four  calibre  tangent  ogive  body  in  transonic  flow  (1^,  =  0.975). 
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This  result  is  taken  from  reference  [13].  The  correction  introduced  into 
the  calculated  surface  pressure  distribution  is  presented  in  Figure  21. 


We  readily  conclude  from  this  result  (i.e.,  Figure  21),  that  provided 
the  boundary  layer  is  thin,  we  can  ignore  the  interaction  between  the 
boundary  layer  growth  and  the  external  potential  flow.  This  result  may 
seem  surprising  in  a  transonic  flow  where  the  nonlinearity  generally 
produces  rapid  changes  in  the  flow  with  small  changes  in  boundary 
conditions.  More  discussions  can  be  found  in  reference  [13]. 


2.  Comments  on  Plume  Induced  Separation 

The  potential  flow  theories  discussed  in  Section  II  cannot 
describe  the  interaction  between  an  exhaust  jet  plume  and  the  body  flow 
field  because  the  details  of  the  interaction  are  intimately  controlled 
by  viscous  forces  in  the  flow.  In  this  section,  some  techniques  for 
describing  the  viscous  effects  are  discussed. 

Five  major  components  of  the  interaction  may  be  isolated.  Each 
will  be  considered,  in  term,  in  the  following  paragraphs.  It  should 
be  emphasized  that  the  applicability  of  the  interaction  studies  presented 
are  restricted  to  supersonic  conditions. 


a .  The  Approaching  Boundary  Layer 

As  a  necessary  initial  condition  for  the  shear-layer  cal¬ 
culation,  the  boundary- layer  development  on  the  body  must  be  known.  It 
is  in  this  way  that  changes  in  the  Reynolds  number  (and  boundary-layer 
transition)  influence  the  interaction  between  the  plurue  and  the  body 
flow  field.  Once  the  flow  is  essentially  turbulent,  however,  further 
increases  in  Reynolds  number  will  be  of  little  significance. 

The  potential  flow  theory  provides  the  basic  pressure  distribution 
from  which  to  calculate  the  boundary- layer  development.  Since  the 
current  theoretical  treatment  is  only  applicable  if  the  flow  is  super¬ 
sonic,  it  is  adequate  to  assume  that  the  presence  of  the  plume  does  not 
influence  the  body  pressure  field  ahead  of  the  interaction.  In  addition, 
it  was  shown  [13]  that  the  distortion  of  the  flow  field  resulting  from 
the  boundary- layer  development  is  also  negligible  for  the  present 
purpose.  Hence,  the  boundary- layer  development  may  be  estimated  directly 
from  the  pressure  field  provided  by  the  theory  of  reference  [ 2] . 

Applying  the  Mangler-Stepanov  transformation  [14,  15]  and  the 
Stewartson- Illingworth  transformation  [16,  17]  to  the  laminar  boundary- 
layer  theory  of  Thwaites  [18]  yields  the  following  result  for  the 
momentum  thickness : 
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Figure  21.  Calculated  Correction  to  Surface  Pressure  Distribution 
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then  the  form  factor  and  skin  friction  coefficient  are  related  to  the 
function  r  by  the  empirical  relations  given  by  Thwaites.  Some  details 
are  also  given  in  reference  [19]. 

Transition  of  the  laminar  boundary  layer  was  considered  to  occur 
either  at  some  prescribed  location  (e.g.,  as  dictated  by  experimental 
evidence)  or  when  the  Reynolds  number, 


Uro&* 


900, 


is  reached  locally  in  the  flow.  This  latter  criterion  was  taken  from  the 
expc  i  ''  ncal  data  presented  in  reference  [20].  If  boundary- layer 
sepLj  'icun  was  predicted  before  the  satisfaction  of  the  above  Reynolds 
number  criterion,  then  the  separation  point  was  used  as  the  point  of 
transition.  Note,  that  separation  is  taken  to  be  when  V  =  -0.09. 

To  complete  the  boundary- layer  calculation,  the  turbulent  boundary- 
layer  theory  of  Nash  [21]  was  used.  This  integral  theory,  based  upon 
the  momentum  integral  equation, 

I  A  (  IT  2  R  \  _  T  .  dUe 

R  dx  (Pe  u  e  R  /  "  TW  •  eue  ° 

is  applicable  to  compressible,  axisymmetric  flow.  The  shear  stress  inte¬ 
gral  is  taken  to  satisfy  a  differential  equation  which  (at  least 
heuristically )  is  representative  of  turbulence  phenomena  (and  is  somewhat 
reminiscent  of  the  equation  of  Burgers  [22]). 

Finally,  if  the  boundary  layer  encounters  a  discontinuity  in 
pressure  distribution  (weak  shock-wave  or  expansion  fan)  then  it  is 
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assumed  that  the  boundary- layer  momentum  thickness  suffers  a  discontinuity 
given  by 


if  L  = 


The  applicability  of  this  equation  to  a  sudden  expansion  (e.g., 
supersonic  flow  over  a  boat  tail  configuration)  is  open  to  speculation. 
Thus,  it  is  shown  in  Figure  22  that  the  boundary-layer  development  on 
such  a  boat  tail  is  extremely  sensitive  to  changes  in  the  initial 
conditions  after  the  expansion.  The  result  quoted  by  White  [23],  being 
of  the  same  form  as  the  above  equation,  is  not  sufficiently  different 
to  give  more  satisfactory  results.  In  addition,  according  to  the  data 
of  Rubin  [24],  there  can  be  a  significant  acceleration  of  the  flow 
ahead  of  the  boat  tail.  The  assumption  of  a  discontinuous  velocity 
change  at  the  boat  tail  is  not  adequate  and  a  more  detailed  theory  for 
this  region  is  required. 


b .  Turbulent  Boundary-Layer  Separation 

In  the  case  of  a  rapid  separation  of  a  turbulent  boundary 
layer  in  supersonic  flow,  the  appr^*  h  ate  theory  of  Mager  [25]  may  be 
used.  Although  designed  for  two-dL  nsional  flow,  the  result  shown  in 
Figure  23  indicates  that  the  pressure  rise  to  separation  in  axially 
symmetric  flow  may  also  be  estimated  by  this  method.  The  experimental 
data  shown  in  Figure  23  is  taken  from  Kuehn  [26]. 

The  data  from  reference  [26]  also  shows  that  above  a  Reynolds 

number  of  R^  ~  10“*,  the  pressure  rise  to  separation  is  little  influenced 

by  increase  in  Reynolds  number.  Indeed,  the  data  in  reference  [27] 

suggests  a  variation  like  Re~l/^0  (which  is  considerably  smaller  than 

the  Rc“^®  variation  predicted  by  Ray  [  28]  X  For  the  present  study,  we 
neglect  the  influence  of  Reynolds  number  on  the  pressure  rise  to 
separation.  Reynolds  number  then  only  enters  the  calculation  as  it 
influences  the  boundary- layer  momentum  thickness  at  the  start  of  the 
interaction.  It  is  considered  that  this  neglect  of  Reynolds  number  on 
the  pressure  rise  to  separation  is  of  higher  order  than  the  accuracy  of 
the  basic  theory  (Figure  23)  for  separation  prediction. 

It  is  also  shown  in  reference  [26]  that  the  sensitivity  of  a 
boundary  layer  to  separation  is  a  function  of  the  velocity  profile.  The 
further  removed  the  profile  is  from  its  equilibrium  value,  the  more 
easily  will  the  boundary  layer  separate.  This  point  will  be  of 
significance  in  the  later  discussion  on  the  flow  past  a  body  with  a 
boat  tail. 


The  pressure  rise  to  separation  also  depends  on  the  downstream 
conditions  (Figure  23)  .  Although  this  result  is  contrary  to  the  usually 
accepted  situation  in  supersonic  separation  (as  embodied  in  the  results 
of  reference  [29]),  it  is  physically  reasonable  to  expect  some  small 
dependence  on  downstream  conditions  resulting  from  the  reversed  flow  in 
the  separation  region.  This  point  is  of  significance  in  the  plume 
induced  separation  problem  where  the  downstream  conditions  (plume  shape 
and  shear  layer  confluence)  depend  upon  the  jet  pressure  ratio.  It  is 
hoped  that  future  developments  can  go  towards  evaluating  the  magnitude 
of  this  "upstream  influence." 


The  equations  of  Paragraph  1  determine  the  boundary- layer  momentum 
thickness  at  the  start  of  the  interaction.  Following  the  suggestion  of 
McDonald  [30],  this  momentum  thickness  may  be  related  to  that  at  the 
start  of  the  shear  layer  by  the  relation: 


v-l  2 

i+y  m2 


[l  -  0.605  Cp  +  11  .725  Cp2  J 


where 


Cp 


1-Mn 


2 

1 


and  Hi  is  the  Mach  number  ahead  of  the  shock,  M2  the  Mach  number  behind 

the  shocks  and  Q  the  boundary-layer  momentum  thickness  at  the  start  of 
the  interaction. 


c  .  Free  Shear  Layer 

It  is  assumed  (following  Kirk  [31])  that  the  free  shear 
layer  leaving  the  body  surface  at  the  separation  point  can  be  replaced  by 
some  equivalent  asymptotic  shear  layer;  the  correct  asymptotic  layer  being 
selected  by  matching  momentum  thicknesses  between  the  shear  layer  and 
the  separated  boundary  layer  at  the  separation  point.  Clearly,  as  the 
free  shear  layer  is  reduced  in  length  (compared  to  its  width)  the 
utility  of  the  equivalent  shear  layer  concept  is  reduced.  In  cases 
where  the  flow  is  not  much  different  from  that  in  a  bubble,  the  concept 
is  entirely  inappropriate. 

Since  we  have  included  the  effect  of  the  boundary  layer  on  the 
shear  layer  development,  the  results  are  more  sophisticated  than  those 
obtained  by  Dixon  et  al  .  [32]  or  Schulz  [33]  and  go,  at  least  some  way, 
towards  including  Reynolds  number  effects. 
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The  flow  equations  in  the  asymptotic  shear  layer  were  described  in 
reference  (1)  and  need  not  be  repeated  here. 

The  unknown  constant  (the  boundary- layer  momentum  thickness  at  the 
start  of  the  shear  layer)  is  determined  from  the  boundary- layer  calculation 
outlined  in  Paragraph  a  and  the  relation  across  the  interaction  given  in 
Paragraph  b. 

d.  Jet  Plume 


Initial  calculations  have  been  performed  using  the  plume 
geometry  developed  empirically  by  Herron  [34].  Further  details  were 
given  in  reference  [  1]  .  At  the  same  time,  a  more  complete  calculation 
using  the  method  of  characteristics  program  developed  by  Prozan  [35]  is 
being  undertaken. 


e.  Conditions  at  Confluence 


We  follow  Korst  and  assume  that  the  final  static  pressure 
after  the  recompression  is  equal  to  the  stagnation  pressure  on  the  dividing 
streamline  before  the  confluence.  Figure  24  shows  this  diagramatically . 

If  it  is  assumed  that  the  recompression  is  isentropic,  then 

7/7-1 

pf  -  Pi(1+2T-Ms) 


where  Pp  and  P^  are  the  static  pressures  before  and  after  the  recom¬ 
pression,  and  Ms  is  the  Mach  number  on  the  dividing  streamline. 

For  the  isoenergic  flow,  we  have: 


with 


Me 


=  (u/Ue>y[l-C^  (u/Ue)2] 


-  7-1  m2 


+  ZlI  M: 


and  u/Ue  is  the  velocity  distribution  in  the  approaching  shear-layer. 
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The  solution  of  the  confluence  region  is  closed  by  means  of  the 
relations  between  pressure  change  and  turning  angle  in  linearized  super- 
sonic  flow.  The  solution  is  valid  only  when  the  flow  downstream  of  the 
surface  separation  is  supersonic  because  this  supersonic  turning  angle 
is  assumed. 


f.  Calculated  Results 


A  typical  missile  configuration  (4-calibre  tangent  ogive 
body)  was  used  to  obtain  data  on  plume  effects  (Figure  25). 

Initial  calculations  were  performed  at  a  free  stream  Mach  number  2 
on  the  body  shown  in  Figure  25  without  the  boat  tail.  The  extent  of 
separation  as  a  function  of  jet  pressure  ratio  is  shown  in  Figure  26. 

The  calculations  were  repeated  at  a  free  stream  Mach  number  1.1. 

Again  the  body  was  taken  without  the  boat  tail.  For  low  jet  pressure 
ratio  (Pj  ~  10  P^),  it  is  seen  in  Figure  27  that  no  flow  separation  is 

present.  At  significantly  higher  pressure  ratios,  separation  is 
produced . 

Figure  28  presents  two  important  points.  First,  it  is  shown  that 
variation  in  the  boundary- layer  transition  point  has  very  little  influence 
on  the  location  of  the  plume  induced  separation  point.  Only  when  the 
transition  occurs  in  the  immediate  vicinity  of  the  separation  point  will 
the  flow  exhibit  a  significant  variation  with  the  actual  location  of 
transition.  Secondly,  Figure  28  shows  that  the  influence  of  the  plume 
has  provoked  separation  of  the  boundary  layei  ahead  of  the  point  of 
separation  predicted  by  the  boundary  layer  calculation  in  the  absence 
of  the  plume. 
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Section  IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  fairly  complete  understanding  of  the  inviscid  transonic  flow  over 
various  bodies  of  revolution  was  obtained  during  this  study.  Engineering 
methods  for  the  computation  of  the  flow  over  an  ogive-cylindrical  body 
with  small  angle  of  attack,  the  conical  nose  with  and  without  angle  of 
attack  and  the  cone-cylindrical  body  at  zero  angle  of  attack  were 
developed.  The  analytical  results  check  very  well  with  the  available 
experiments.  This  Increases  our  level  of  confidence  in  the  analysis 
discussed  herein. 

Many  improvements  are  still  needed  to  obtain  more  general  solutions. 
For  instance,  application  of  the  existing  theory  to  an  arbitrary  shaped 
smooth  body  (including  the  acceleration  and  deceleration  portions)  is 
very  desirable.  The  improved  solution  around  the  apex  point  and  the 
shoulder  for  a  cone-cylindrical  configuration  should  also  be  studied. 

To  include  the  angle  of  attack,  or  the  effect  due  to  the  cross- flow, 
on  the  local  two-dimensional  method  is  another  subject  needing 
careful  study. 

It  was  found  that  the  viscous  effects  on  an  ogive-cylindrical  body 
are  not  too  significant.  However,  estimation  of  the  viscous  effects 
applicable  to  the  boat  tail  portion  needs  specific  attention.  Our 
preliminary  viscous  analysis  indicates  that  the  shoulder  of  a  boat  tail 
is  very  crucial  to  the  boundary  layer  separation  caused  by  the  exhaust 
jet  plume. 

The  calculation  of  an  entirely  subsonic  or  a  mixed  supersonic/ 
subsonic  "separation  will  form  the  subject  of  further  study. 
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